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ABSTRACT 

We present here X-ray spectra of the HMXB SMC X-1 obtained in an observa- 
tion with the XMM-Newton observatory beginning before echpse and ending near 
the end of echpse. With the Reflection Grating Spectrometers (RGS) on board 
XMM-Newton, we observe emission hnes from hydrogen-hke and hehum-hke ions 
of nitrogen, oxygen, neon, magnesium, and sihcon. Though the resolution of 
the RGS is sufficient to resolve the helium-like n=2— i>l emission into three line 
components, only one of these components, the intercombination line, is detected 
in our data. The lack of flux in the forbidden lines of the helium-like triplets is 
explained by pumping by ultraviolet photons from the BO star and, from this, 
we set an upper limit on the distance of the emitting ions from the star. The 
lack of observable flux in the resonance lines of the helium-like triplets indicate 
a lack of enhancement due to resonance line scattering and, from this, we derive 
a new observational constraint on the distribution of the wind in SMC X-1 in 
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velocity and coordinate space. We find that the sohd angle subtended by the 
volume containing the helium-like ions at the neutron star multiplied by the ve- 
locity dispersion of the helium-like ions must be less than 47rsteradianskms~^. 
This constraint will be satisfied if the helium-like ions are located primarily in 
clumps distributed throughout the wind or in a thin layer along the surface of 
the BO star. 



1. Introduction 

In isolated early-type stars winds are driven as ultraviolet photons from the stellar 
surface impart their outward momentum to the wind in resonance line transitions. In a 
high-mass X-ray binary (HMXB), the wind is ionized by X-rays from the compact object. If 
the X-radiation is intense enough, the resulting ions will not ha ve transitions in the ultraviolet 



and this greatly affects the dynamics of the wind (see, e.g., lBlondinlll994l . and references 
therein). The SMC X-l/Sk 160 system, which consists of a 0.71 second X-ray pulsar and a 
BOI star together in a 3.9-day orbit, is the most X-ray luminous known HMXB and, therefore, 
presumably, an extreme example of wind disruption by X-ray ionization. 

The behavior of the wind of an early-type star under the influence of ionization from an 
X-ray emitting companion has been the subject of many theoretical studies. For SMC X-1 , 



the most relevant and detailed such study is, arguably, the one by iBlondin fc Wool (119951 ) 
which included a numerical hydrodynamic simulation. The features that appeared in that 
simulation included a regular, UV-driven wind on the X-ray shadowed side of the star, 
a thermal wind on the X-ray illuminated side, transition regions between the two types of 
winds and dense, finger-like structures in the equatorial plane. However, even this simulation 
includes significant approximations: the gravity of the compact object is not included and 
X-ray photoionization and its dynamical effects are treated in a very approximate way. 
Because of the complexity of these systems and the difficulty of accounting for all of the 
relevant physics, observations are critical to characterizing the behavior of HMXB winds. 

Though previous X-ray observations of SMC X-1 have had spectral resolving powers 
(= A/AA = E/AE) of ~50 or less. X-ray spectroscopy has already revealed much about the 
wind. X-ray spectroscopic observations of this resolving power can, through measurements of 
line fluxes, indicate the quantity and ionization level of X-ray emitting material. In observa- 
tions of SMC X-1 with ASCA, emission lines were not detected and upper limits were set on 
the quantity of material of intermediate ionization level in the wind of SM C X-1, excluding 



the presence of finger-li ke structures as they appeared in the simulation of iBlondin fc Woo 



(IWojdowski et al.ll2000l ). Further observations of SMC X-1 with the Advanced Camera for 
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Imaging Spectroscopy (ACIS) on-board the Chandra X-ray Observatory, however, have de- 
tected emission hnes a nd, therefore, the presence of material of intermediate ionization in 
the wind of SMC X-1 ( IVrtilek et aLll200ll ). These results indicate that overdense r egions do 
exist in the wind of SMC X-1 and that while the simulations of Blondin &: Wool ma.v not 
be accurate in detail, the general form of the structure predicted by iBlondin fc Wool may, 
in fact, be present. To develop further constraints on models of the wind in SMC X-1, it 
is desirable to have constraints in addition to the quantity and ionization level of the X-ray 
emitting plasma. 

High-resolution spectroscopic observations with resolving powers >200 have the poten- 
tial to provide useful constraints on the kinematics of the high-ionization wind in HMXBs 
through direct measurements of Doppler line shifts and broadening. In addition, the n = 
2^1 triplets of He-like ions are resolved at high resolution and measurements of the fluxes 
of the individual lines of these triplets provide a constraint on the structure and kinematics 
of the X-ray emitting material. Of the three components of the helium-like triplet, only 
one — the resonance line — has a large oscillator strength and can be enhanced by reso- 
nant line scattering. Because resonant line scattering saturates and because this saturation 
depends on the distribution of the scatterers in physical and velocity space, constraints on 
the distribution of the wii id in physical and velo city space can be derived from the relative 



fluxes of these three hnes (IWojdowski et al 



20031). The flux ratios of the heli um-like triplets 



are also affected by ultraviolet radiation ( IBlumenthal. Drake. &: Tuckerll 19721 ) and, therefore, 
the flux ratios measured from an HMXB constrain the distance from the photosphere of the 
high-mass star to the emission region. 

Because of the location of SMC X-1 outside of the plane of the Galaxy, the column 
density of interstellar material to it is small and, unlike the HMXBs in the Galaxy, it can 
be observed in the wavelength range 15-35 A. The n = 2 — 1 transitions of the hydrogen- 
like and helium-like ions of oxygen and nitrogen as well as hydrogen-like carbon lie in this 
range. Therefore, unlike with the Galactic HMXBs, it is possible to derive observational 
constraints on the density and velocity distributions of the regions of moderate ionization in 
SMC X-1 where these ions exist. Because of the high spectral resolution (FWHM of ~65 mA, 
corresponding to resolving powers in the range 230-540) and large effective area (total of 
~80cm^) of the Reflection Grating Spectrometer on XMM-Newton in this band, the RGS 
on XMM-Newton is particularly well-suited to spectroscopic study there. We present here 
observations of SMC X-1 with XMM-Newton beginning before an eclipse and ending near 
the end of that eclipse. We present an analysis in which we focus on the data from the RGS 
and derive constraints on the space and velocity distribution of material in the wind. In 
^we describe our observations and an empirical analysis of the line emission in which we 
measure line fluxes, shifts, and widths. In ^we derive constraints on the distribution of the 
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wind in physical and velocity space from the line fluxes we measure. In §U we discuss the 
implications of our results for models of HMXB winds. 



2. Observations & Empirical Analysis 

SMC X-1 was observed from on 2001-05-31 from approximately 02:14 to 18:50 UT with 
XMM-Newton. The actual times at which the first and last observational data were received 
with the various instruments differ from those times by as much as an hour. However, data 
were collected with each instrument during at least 90% of the time between the nominal 
beginning and end of the observations. In Figure [T] we show plots of count rates from 
the source region and a background region on the MOSl detector. In the beginning, the 
count rate in the source region decreases through X-ray eclipse ingress. Toward the end of 
the observation, the count rates for both the source and the background regions increase 
while the difference stays approximately constant indicating an increase in the non-X-ray 
background. We divided the data into the three time intervals noted on the figure. On the 
time scale of the figure, the times which divide the intervals are 10609 s and 42289 s. The 
intervals, labeled "ingress", "center", and "end" were chosen to divide the eclipse ingress 
from the time of complete X-ray eclipse and the time of high background from the rest 
of the observation. We then extracted spectra from each of the five XMM-Newton X-ray 
instruments for each of the three intervals and from the interval consisting of the center and 
end intervals which we label "eclipse" . 

We focus here on the spectral data from the RGS. However, in Figure [2] we plot the 
spectrum from MOS 1 fo r the eclipse interva l. This spectrum is similar to those spectra 



of SMC X-1 obtained by IVrtilek et all (120011 ) with Chandra/ ACIS. In Figure E] we show 



our RGS spectrum for the eclipse interval. Line emission from hydrogen-like and helium- 
like ions of several elements is apparent. The spectral resolution of the XMM-Newton RGS 
is such that the n = 2 ^ 1 transitions of helium-like ions are resolved into a triplets 
consisting of, in order of wavelength, the resonance {ls2p^Pi aground), intercombination 
(ls2p^P2,i —aground), and forbidden (ls2s —aground) lines. However, in our RGS spectra, 
only the intercombination lines are apparent. Before attempting to explain this fact or to use 
it to derive constraints on the properties of the SMC X-1 system, we first make a quantitative 
characterization of the line emission. 

In order to measure the fluxes, widths and shifts of these lines, we fit the wavelength 
channels within 0.5 A of the rest wavelengths of each of the n = 2 —>■ 1 emission line 
complexes. We fit the spectrum in the vicinity of each line complex using a power law 
continuum and Gaussian lines. We express the width (sigma) and centroid of the Gaussians 
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Fig. 1. — Count rates from the MOSl detector. We show the count rate from a source 
region in black and from a background in red. An increase at late time in both count rates 
is indicative of an increase in the non- X-ray background. Each data point represents the 
average count rate over 128 s. Time zero on the abscissa is 2001-05-31 02:22:16.9357 UT. 
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Fig. 2. — MOS 1 Eclipse spectrum of SMC X-1. Labels indicate the energies of spectral 
features which may or may not be detected in this spectrum. The labels denote the energies 
of lines except where the ion is followed by "e" indicating an edge. The line energies indicated 
are n = 2 ^ 1 transitions except for VIII at 0.77 keV, Ne X at 1.21 keV, and Mg XII at 
1.74 keV which are n = 3 — 1 transitions. 
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in velocity coordinates and fix the width and shifts of the hnes in a single complex to be the 
same, i.e. : 

^ Xk{vr/c)^/^^ \ 2Xl{o^^rlcy J 

where Fx is the photon (not energy) flux, the A^s are the rest wavelengths of the lines, and 
a, fr, o"^,r, and the /^s are fit parameters. We assume that lines appear in emission and 
therefore constrain the fluxes of all lines to be greater than or equal to zero. 

The inter combination line is actually an unresolved doublet. Therefore, we fit the ingress 
and eclipse spectra in the vicinity of the helium-like triplets using the model described by 
equation [T] with four Gaussians: one for each of the resonance and forbidden lines and one 
for each of the two components of the intercombination line. However, we fix the ratio of the 
fluxes of the two compon ents of the intercombination line to be equal to that predicted by 



Kinkhabwala et al.l (120031 ) for recombination using temperatures given in Table [31 Further- 
more, for reasons we describe later, we formulate our model not in terms of the individual 
fluxes of the three lines of the triplets but in terms of the total triplet line flux / and the 
ratios of the line fluxes G and R where 

(2) 

and 

R^j, (3) 

where the subscripts "r" , "i" , and "f" indicate, respectively, the resonance, intercombination, 
and forbidden lines. In the fits, we restrict the value of a^^r to be less than 1500 km s~^, Vr 
to be between —1000 and 1000 km s~^ and we restrict J, G, i?, and K to be greater than or 
equal to zero. In Figures H] and [5], we show our best fits to these triplets and in Table [1] we 
show the best-fit values and confidence intervals for these line fluxes, shifts and widths. The 
fact that we detect only the intercombination lines of the helium-like triplets is reflected in 
the fact that we obtain only lower limits on G and only upper limits on R. 

The hydrogen-like n = 2 1 transitions consist only of an unresolved doublet and we 
fit these doublets using two Gaussians with the ratio of the fluxes of those two lines fixed 
at 2:1, proportional to the statistical weights of the upper levels and the ratio expected for 



^This awkward model for the continuum flux is due to the fact that the power law spectral model is 
defined in energy space {Fe = K{E /keV)"). Our fits to the contimmm flux arc not a focus of this work nor 
do we make any use of them here so we proceed with this continuum model. 
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Fig. 4. — Helium-like triplets for the "ingress" interval. 
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Fig. 5. — Helium-like triplets for the "eclipse" interval. 



Table 1. Measured Helium-like Triplet Parameters 



power-law 



Element 


ja 


G 


R 


(T^,r(knis ^) 


i'r(kms ^) 


K 


a 




ingress 


N 



Ne 

Mg 

Si 


^•"-1.5 

-'-•'±-1.4 
9+6 

22+^^ 


>4 
>3 


<0.3 


<800 
5001^°° 


-3001^°° 
-500ir 


170±30 
2051^6 

280+?° 


u.y_2 

-'-•-'--0.3 

-1.3 ±0.5 


eclipse 


N 



Ne 

Mg 

Si 


0.22±0.18 
1.1±0.4 

9 9+0.9 
^■^-0.7 

2 4+2-^ 


>1.6 

>9 
>2 


<0.63 

'-'•^-0.2 

<0.25 


<500 
<300 
<600 
<1300 


_4QQ + 500 

gQ+130 

_4QQ + 500 

-6001^° 


<1.3 
<0.0002 
6.5±0.7 

Q Q+1863.5 

9 r+1863.9 
^■^-0.8 


91 c+o.e 

^-■^■^-13.6 

-6lL 

01- 
01- 



Note. — As a confidence limit on a^^r, v^, or a, oo indicates that parameter boundaries 
are within the confidence limits. The entry " • • • " indicates that the data do not constrain 
this quantity. 

*Total triplet line flux in units of 10~^ photons cm~^s~i. Not corrected for interstellar 
photoelectric absorption. 
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recombination. Furthermore, we restrict to be less than 2500 km s~^ and Vr to be between 
—2500 and 2500 km s~^. In Table [2], we show the fit parameters for the hydrogen- like lines. 
Again, / indicates the total photon flux for the doublet. 

3. Line Fluxes as Measures of the Distribution of Scatters in Physical and 



The emission lines from high-mass X-ray binaries result from recombination and reso- 
nant line scattering in wind of the high-mass star. The line emission from these processes 
depends on the distribution of the wind material in coordinate space and, in the case of 
resonant line scattering, in velocity space. Therefore, the line fluxes we measure provide a 
constraint on these distributions of the wind material. However, the equations that describe 
these distributions are complex and rather than attempt to derive constraints directly from 
the line fluxes, we believe it will be more illuminating to create a simple model described 
by a few parameters and then derive constraints on these parameters. In §3.11 we derive 
equations for line luminosities due to recombination and resonant line scattering. In §3.2[ 
we describe a modification to those luminosities for the case of the helium-like triplets due 
to ultraviolet radiation of the high- mass star. In §3.31 we describe the simple models and our 
constraints on their parameters. 



In a plasma in photoionization equilibrium, photons are absorbed from the continuum in 
photoionizing transitions and line photons are emitted as electrons and ions recombine. If the 
plasma is optically thin in the photoionizing continuum, the recombination line luminosity 
may be expressed as 



where h is Planck's constant, uq is the frequency of the line, e and me are the electronic charge 
and mass, nion is the density of the ion that emits the line, rj is fraction of recombinations 
to the ion which result in emission of the line, Lj, is the specific luminosity of the compact 
X-ray source, i/th is the ionization threshold frequency of the emitting ion in its ground state, 
and df/du is its continuum oscillator strength, and s and fl are the distance from and solid 
angle subtended at the compact X-ray source. 

In any plasma, line photons are absorbed and reemitted in resonance line transitions. 



Velocity Space 



3.1. Recombination and Resonant Line Scattering 
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Table 2. Measured Hydrogen-like Line parameters 



power-law 



z 


/(10-V^cm-i) 


cr^,r (km s ^) 


Vr(kms ^) 


norm 


a 




ingress 


6 

7 

8 

10 

12 

14 


'-'•'-'-0.5 

4li 
7±5 
1.4l?i 
441^^ 


<1200 

^'-"-'-900 

25001^00 


-2300ir 
20001^00 


<1900 
17411L 

17lli^ 
170±30 
290001™ 


oir 

-611, 

o.ol~ 

-0.4lo°°3 

-8.611:^ 


eclipse 


6 
7 

8 

10 
12 
14 


0.3±0.3 
0.7±0.2 
1.7±0.3 
1.0±0.6 

<2.0 


13001^0 
<400 

<400 


-15001^°° 
0±200 

2QQ+900 


<8 

9±3 
8.4±1.9 
7.5±1.4 

17+1892 
q q+2632.6 

'J-'J-i.e 


olr 

oir 



Note. — The entry " • • • " indicates that the data do not constrain this 
quantity. 
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The net effect of this is only to redirect hne photons. However, depending the geometry 
of the plasma and radiation source, resonant line scattering may result in absorption or 
emission lines being observed. For the geometry of interest in this work — an occulted 
compact radiation source with radiation scattered into the observer's line of sight by plasma 
— resonant scattering results in emission lines and the line luminosity due to scattering may 
be expressed as 

L,„(47r)-^ /"/(I - e-^^'''"^)diydn (5) 



^scat 



where r is the line optical depth through the plasma along a line from the compact radiation 
source. The line optical depth is given by 

Tre^ r 

T = f / nion(f){iy)ds (6) 

rrieC J 

where / is the line oscillator strength and the function (p describes the line profile. Neglecting 
natural line broadening and considering only Doppler line broadening, we have 

cPiu) ^ -g{c{u/uo - 1)) (7) 

where g is the local distribution of the ions' velocity component in the direction along the 
line of sight from the compact X-ray source and is normalized such that 

g{v)dv = 1. (8) 

However, we assume that ion velocities are very small compared to c so that g{v) is non-zero 
only for \v\ <^ c. With this definition, 

(f){iy)du = 1. (9) 



As we will ultimately be interested in calculating observed line fluxes the spatial integrals 
in equations H] and O should be restricted to the volume that is both illuminated by the 
compact radiation source and observable (i.e., not occulted by the star). In principle, the 
radiation from the compact source that reaches the observable plasma may be depleted at 
the line frequency by line scattering in the occulted region, reducing the observable line 
luminosity relative to that described by equation [51 However, as scattering does not destroy 
photons but merely redirects them, for most plasma geometries, the total number of line 
photons reaching the observable region will not be significantly affected by scattering in the 
occulted region. Only for a specialized plasma geometry: a long, narrow distribution of 
plasma, oriented along a line of sight from the compact radiation source will equation [5] be 
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significantly inaccurate. Furthermore, if the part of the wind where the hne emission occurs 
is moving supersonically (as HMXB winds are known to do) the hne frequencies in different 
parts of the wind will be shifted apart so that resonant line scattering in one part of the 
wind will not be affected by resonant line scattering in another part. 

In equation HJ we have t] inside the spatial integral, allowing for the possibility that 
the recombination efficiencies may depend on the local conditions. In fact, for helium-like 
lines, the forbidden line may be "pumped" into the intercombination line at high densities 
or high ultraviolet radiation intensities and we define rj so that the effects of this pumping 
are absorbed into it. Near a hot star such as Sk 160, the BO companion of SMC X-1, 
the ultraviolet radiation intensity is high. Therefore, the line luminosities due to resonant 
scattering depend on the distance from the star. In the next section, we describe this in 
more detail. 



3.2. Helium-like Triplets 

As can be seen from the equations above, the observed flux of any one line depends in 
a complicated way on the distribution of the emitting ions in both coordinate and velocity 
space. Therefore, the constraint that can be derived from a single line is quite complex. 
However, because the resonance line of each helium-like triplet has a significant oscillator 
strength while the other two lines do not, only the resonance line is affected by resonant 
scattering. Therefore, only the flux of the resonance line depends on the velocity distribution 
of the emitting ions and by measuring the flux of all three lines of the triplet, it is possible to 
derive more specific constraints on the distribution of the emitting ions in space and velocity 
than can be obtained from the flux of any one line. 

As previously mentioned, the recombination efficiencies of the intercombination and 
forbidden lines are not constant but depend on the local conditions. This is due to the 
fact that the upper level of the forbidden line, the ^5*1 level, is metastable and an ion 
in that state may be excited by absorption of an ultraviolet photon or collision with an 
electron to any of the levels. Two of these levels {^P2 and ^Pi) make up the upper 
levels of the intercombination line. Therefore, at high densities or high ultraviolet radiation 
intensities, the effective recombination efficiency of the intercombination line is increased and 
the effective recombination efficiency of the forbidden line is decreased. This effect has been 
used by several authors to estimate the distance of emitting ions from isolated hot stars. 
However, several of the calculations have neglected the fact that once an ion is excited from 
the ^5*1 state to a ^P state, it may decay back to the ^5*1 state by emission of an ultraviolet 
photon rather than decay to ground by emission of an intercombination line photon. We 
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write the recombination line efficiencies including this effect as 

r/i = r^^ + r/Kl-fef) (10) 

Tjf = 7][hi (11) 

where ri[ and ri[ are the recombination efficiencies without ^Si — > excitation and 6f is 

the branching ratio for radiative decay of the ^Si state to groun d. In order to der i ve an 



expression for this branching ratio, we solve the rate equations of iMewe fc Schrijverl (11978 



equations 18 and 19) and neglect inner-shell ionizatio n, all collisions, and in duced radiative 



decay. To write our expression, we use the notation of lPorquet et al.l (1200 ll ) and denote the 



^S'l level by m (for "metastable"), the ground (Is^ ^Sq) by g and denote the level ^Pk by pk 
where k may be 0, 1, or 2 and get 

h, = — (12) 
4 _l_ T 



A -i- A 

k=Q ^Pk9 ^ ^Pkm 



where 



where Wmp^. is the excitation rate from the to ^PjU. 
For photoexcitation, the w values are given by 

w =^f ^ (14) 

where Ju„,p^ is the mean intensity of radiation (i.e., the intensity averaged over solid angle) at 
the frequency Ump^ and Umpk and fmp^ are, respectively, the frequency and oscillator strength 
of the transition m pk- The frequencies i^mpk for elements from carbon to iron range from 
the near to the extreme ultraviolet. Because the quantity rj^ + rjf and the value of ?7r do not 
depend on the ultraviolet pumping, neither does the value of the G ratio but the value of 
the R ratio does. 



^Our expression for the f o rbidd en line e fficiency is cons i stent with expressions for forbidden fine intensities 
given bv lMewe fc Schriiverl (jl978r ) and bv lPorguet et alj (j200ir ). However, our expression for the intercom- 
bination line is not consistent with those works. As they are written, it is difficult to compare our line 
efficiencies with those line intensities. Therefore, we note that our sol ution of the rate equa tions implies 

— J2k=i ^gpki^ ~ BRp^mBRmg) + Cgm{^ — BRmg) for cquatiou 5 of lPorquet et al.l (|2001l ). We are en- 
couraged in using our expression for the intercombination line because, with our expression, the sum of the 
intercombination and forbidden lines, as expected, does not depend on the magnitude of the '^S'l -^^ P 
excitation. In any case, the final results for both of those works are based not on those expressions but on 
solutions of the full rate matrix and would not be affected by any error in the expressions. 
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3.3. Simple Model 



From the observed fluxes of tlie tliree lines of a given helium-like triplet, it is possible to 
derive constraints on the distribution of the helium-like ions in coordinate and velocity space. 
However, because of the complexity of the equations describing the line luminosities, the 
constraints that can be derived are also complex. Therefore, we compute line fluxes for two 
simple parameterized models for the distribution of material. Then we derive constraints on 
those parameters from the observed line fluxes. We do not actually expect the distribution of 
material in the wind of SMC X-1 to be described by these simple models in detail. However, 
from these model parameter constraints it is possible to make approximate inferences about 
the distribution of the wind of SMC X-1 in coordinate and velocity space. 

In our models, each ion exists exclusively within a solid angle subtended at the neutron 
star, and is distributed such that the column density along lines of sight from the neutron 
star has the single value N within the solid angle Q. The distribution of ion velocities along 
lines of sight from the neutron star in our model is uniform over some range with width 
Av (which need not be centered on zero). Also, to simplify the calculation of ultraviolet 
pumping, we take the entire plasma to be at a single distance r from the surface of the 
companion star which we take to be a sphere emitting as a blackbody with a temperature of 
=30,000 K, the approximate effective temperature for a star of type BO such as Sk 160. 

For this model, the UV radiation density is 



J, = W{r)B,{T^) 



(15) 



where is the Planck function. 
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and is the stellar radius. The luminosity of a hue due to recombination is then 




(17) 



and the line luminosity due to scattering is 



L,,,, = u,{i7T)-'^nAv{i 




(18) 



or, equivalently. 



1 




N/Av 



(19) 
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It can be seen from the above expressions that it is not necessary to specify all of the 
model parameters in order to determine the line luminosities. It is sufficient to specify only 
the parameter combinations QN and N/Av rather than the values of the three parameters 
contained in those two combinations. If we replace the three parameters with the two 
parameter combinations then the fluxes of the lines depend linearly on QN, the sum of the 
luminosities of the intercombination and forbidden lines depends only on QN, the G ratio 
depends only on N/Av, and the R ratio depends only on r. 

In principle, it should be possible to derive constraints on the model parameters from 
the constraints on the line fluxes in Table [H However, for the line emission mechanisms we 
consider here, the G a nd R ratios can only take on values within a finite range (e.g., see 



Wojdowski et al.ll2003l Figure 8 for the range of the G ratio for Si XIII). In at least one case, 
our best fit value of G falls outside of this range. Therefore, in order to derive meaningful 
constraints on the physical model parameters, we redo the fits with the G and R values 
constrained to be within the range allowed by scattering and recombination as described 
above. In all cases, we are able to obtain good fits with these constraints imposed. 

In order to relate our model parameters to line fluxes and derive constraints on our 
model parameters, we use the following data. We get recombination efficiencies (without 
P pumping) by dividing effective line recombination rates from iKinkhabwala et al. 



(I2OO3I ) by total recombination rates from IVerner et al.l (Il996l ). We have written equations 



for line luminosities in terms of the specific luminosity of the X-ray source. Of course, we 
measure fluxes and any luminosities we derive are subject to uncertainty in the distance to 
the object. However, if we write equations H] and [5] in terms of fluxes, factors of distance 
cancel. Therefore, in our calculations we use fluxes and our results do not depend on the 
distance to SMC X-1. For the specific luminosity of the SMC X-1 neutron star we use 



where JF is the specific photon flux per unit photon energy interval which is related to L^, 
by Li, = {4:TxdP)^^h'^vJ^ and e is photon energy with K = 0.0267 photons cm~^ ^^Y"^ 



a = 1 from a fit to the uneclipsed ASCA spectrum of SMC X-1 (I Wojdowski et al.ll2000l ). 



The value of 1 we use for a is an approximation to the best-fit value of 0.94±0.02. For the 
continuum oscillator strength of the helium-like ions above the ionization threshold we use 



(c.f., IWojdowski et al.ll2003l ). For all elements, we take the oscillator strength of the res- 
onance line to be 0.7, which is accurate to 5%. For line frequencies (wavelengths) and 
spontaneous decay rates we use values from the Astrophysical Plasma Emission Database 
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(APED. ISmith et al.ll200ll ). We take the oscillator strengths fmpo to be zero and derive the 
oscillator strengths fmpi and fmp2 from the spontaneous deexcitation rates of APED. We 
give the results of our model parameter constraints in Table [31 



4. Summary /Discussion 

We have observed SMC X-1 in eclipse with XMM-Newton and, with the RGS, resolved 
the helium-like triplets of nitrogen, oxygen, neon, and magnesium. To our knowledge, this 
is the first time helium-like triplets from SMC X-1 have been resolved and the first time the 
helium-like triplets of nitrogen, oxygen, or neon have been resolved for any HMXB. In all 
cases, only one component of the triplet, the intercombination line, is observably present. The 
lack of observable forbidden line fluxes in the triplet is easily explained by photoexcitation 
pumping of the forbidden line into the intercombination line by the ultraviolet radiation of 
the BO star Sk 160 and allows upper limits to be set on distance of the emitting helium-like 
ions from Sk 160. The absence of observable fluxes in the resonance lines is consistent with 
what we expect from recombination in the photoionized wind. However, the fact that the 
resonance lines are not enhanced by resonant scattering implies a lower limit on the optical 
depth of the wind in the resonance line and, therefore, constraints on the structure and 
kinematics of the wind are also implied. 

We have modeled the helium-like line emission as recombination and scattering emission 
from a region that is partly obscured by the companion star and is described by a single solid 
angle (Q) and column density (A^) and the velocities along lines of sight from the neutron 
star have a boxcar distribution with width Av for each of the ions. We have inferred infer 
the quantity QAN from the sum of the flux of the intercombination and forbidden lines 



Table 3. Simple Model Parameters Derived from Helium-like Line Emission 



element 


T,„c(eV) 


At 




N/Av 
{10l5cm~2(kms-i)-i) 


Q{An)-^N 
(10^5 


n{47r)-iA?) 
(kms-l) 


'Tree 

reference 


N 


3.0 


>4.2 


<57 


> 1.0 


0.22±0.18 


<0.39 


1 


O 


4.0 


>20 




> 6.4 


1.3±0.4 


<0.27 


1 


Ne 


12 


>22 


<6.6 


> 11 


r, r, + 0.9 

^■^-1.0 


<0.27 


2 


Mg 
Si 


21 

75 


>4.2 


<1.4 


> 3.2 

5 1+5.1 


3At\:t 


<1.5 


2 
2 



'^r is the distance from center of the star. 

References. — Q'l iKinkhabwala et al]|2003l: (2) ISako et"al]ll999l 



- 19 - 



which are not affected by resonant hne scattering. We infer the quantity N/ Av from flux of 
the resonance hne relative to the sum of the other two (the inverse of the G ratio). From 
the constraints on these two quantities, we also obtain a constraint on the quantity QAv. 
Because we do not detect the resonance line, we obtain only lower limits on N/ Av and only 
upper limits on QAv. 

As previously mentioned, in our analysis we have assumed that photons escape isotropi- 
cally along lines from the location of their production or first scattering. While we do expect 
recombination emission to be isotropic, resonant line scattering is not isotropic. The angular 
distribution for resonant line scatteri ng is a linear combiii ation of an isotropic distribution 
and an dipole distribution (see, e.g.. IChandrasekharl 119601 . Ch. 19). Furthermore, for large 
optical depths that are comparable to or greater than unity such as we infer, photons may 
undergo several scatterings before escaping. In this case, photons will escape preferentially 
along directions where the optical depth is least. In general, the calculation of photon escape 
direction is also complicated by the fact that photons may travel large distances across the 
plasma before escaping. While it is difficult to include this possibility in analytic calculations, 
it is straightforward using Monte Carlo techniques. However, if the region that emits the 
helium-like line emission has supersonic velocity differentials, as we know that the bulk of the 
wind does, the tra nsfer of resonant line photons is essentially local (the Sobolev approxima- 
tion, ICastorlll97Cll ) and can be approached analytically. For a supersonically expanding wind. 



the optical depth in a direction given by the coordinate x is inversely proportional to dvx/dx. 
If the flow moves along lines from a single point and y is the distance from that point, then 
these derivatives, in directions parallel and perpendicular to the flow, respectively, are dv/dy 
and v/y. While these quantities will not generally be equal, they will generally be of the 
same order over the bulk of the wind and, therefore, if we assume complete redistribution in 
freq uency and dir ection for individual scatterings, then the radiation will escape isotropically 



see 



Castoiill970l ). In fact, individual scattering events do not redistribute frequency and 



direction completely and so even with v/s = dv/ds, radiation does not escape isotropically. 
However the differe nce fron a the case of co mplete redistribution is no greater than 10%-20% 
( ICaroff et al.lll972l . see also lMihalaslll980l ). In light of these facts, we proceed to explore the 
implications of our results. However, it must be kept in mind that the approximations we 
have made for the radiative line transfer may quantitatively affect our conclusions. 

The implication of our results on fundamental wind parameters, such as the mass- 
loss rate and the terminal wind velocity, is complex owing to factors such as the complex 
dependence of the ion fractions on density in photoionization balance. Therefore, interpreting 
these results in terms of fundamental wind parameters is difficult to do without computing 
line fluxes for complete wind models which is beyond the scope of this work. However, the 
upper limits on {Q/47i)Av that we obtain, approximately lkms~^, are quite small compared 
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to 1000 km s~^, the order of the terminal win d velocities of massive st ars and the wind velocity 
in a simulation of the wind of SMC X-1 by iBlondin &: Wool (Il995l ). Therefore, the part of 
the wind that emits the helium-like lines cannot be homogeneously distributed around the 
entire wind. Instead, we believe that our results indicate either that the volume (or volumes) 
containing the helium-like ions of the various elements subtend a small solid angle at the 
neutron star or that the helium-like ions exist primarily in a region where the wind is at a 
small fraction of its terminal velocity or, perhaps, some combination of both. 

One possible explanation of our results is that the helium-like ions exist primarily in 
dense clumps throughout the wind. Dens e clumps in t he wi nd of an HMXB have previously 
been invoked in the case of Vela X-1 by ISako et al.l (119991 ) in order to explain the strong 
fluorescence lines observed from that system. Another possible explanation is that helium- 
like ions exist mainly near the surface of the star where the density is higher, the ionizatio n 
less, and the velocity lower than in the outer parts of the wind (see, e.g.. lLiedahl et al.ll200ll ). 
Furthermore, the part of a volume along the surface of the companion star with a thickness 
signiflcantly less than one stellar radius that is visible during eclipse would subtend a solid 
angle signiflcantly less than 4it from the neutron star. In Figure [6] we illustrate both of 
these possibilities. Deflnitive tests of these hypotheses would require calculations of the line 
emission from detailed wind models and that is beyond the scope of this work. However, we 
expect that our results, including the remarkably small values of flAv are consistent with 
current expectations about the nature of winds in HMXBs. 
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Fig. 6. — In this figure we illustrate to possible distributions of the helium-like ions consistent 
with the helium-like line fluxes we measure. In the first panel we illustrate the possibility 
that the helium-like ions are in clumps distributed throughout the wind. In the second panel 
we illustrate the possibility that the helium-like ions are near the surface of the companion 
star. 
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